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ABSTRACT 



Aims. We search for periodic variations in the radial velocity of the young Herbig star HD 135344B with the aim to determine a 
rotation period. 

Methods. We analyzed 44 high-resolution optical spectra taken over a time range of 15 1 days. The spectra were acquired with FEROS 
at the 2.2m MPG/ESO telescope in La Silla. The stellar parameters of HD 135344B are determined by fitting synthetic spectra to the 
stellar spectrum. In order to obtain radial velocity measurements, the stellar spectra have been cross-correlated with a theoretical 
template computed from determined stellar parameters. 

Results. We report the first direct measurement of the rotation period of a Herbig star from radial- velocity measurements. The 
rotation period is found to be 0.16 d (3.9 hr), which makes HD 135344B a rapid rotator at or close to its break-up velocity. The rapid 
rotation could explain some of the properties of the circumstellar environment of HD 135344B such as the presence of an inner disk 
with properties (composition, inclination), that are significantly different from the outer disk. 
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1. Introduction 



; Herbig Ae/Be /HAeB e ) stars (iHerbid IT960t 
■ iFinkenzeller & Mundtl Il984t IWaters & WaelkensI Il998l) are 
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intermediate-mass (2-10 Mq) pre-main sequence stars. They 
are believed to be higher-mass counterparts of the T Tauri stars 
(TTS) and, therefore, fill the parameter space between TTSs 
and high-mass young stars in addressing the que stion of star 
formation as a function of mass (I AppenzelleJ 1 9941) . 
Stellar rotation is a crucial parameter for the evolution of 
angular momentum, magnetic fields, and accretion processes. It 
is known that HAeBe stars exhibit significantly larger projected 
rotational velocities, vsin/, than TTSs. Typical v sin / values for 
HAeBe stars are in the range of 60 to 225 km s ~^ (iDavis et alJ 
[19831 lFinkenzelIeii[T985l iBoehm & Cataial[T99l . whereas most 
TTSs have vsin; values of about 10 km s~' (e.g. LWeLse^et alj 
I2010I) . This result indicates that mechanisms of angular momen- 
tum d ispersal work much les s efficient in HAeBe stars than in 
TTSs dBoehm & Catalal[l995h . Star-disk locking and subsequent 
rotational braking might even be absent in HAeBe star s. Only 
for low -mass HAeBe stars (M < 2.6 Mq) Boehm & Catalal 
(Il995h found indications for loss of angular momentum due to 
stellar winds. 

However, there is no extensive study to determine the rotation 
periods of HAeBe stars present yet. Only for a few HAeBe 
stars, rotation periods were found so far. Observed vaiiations 
of Ca II K and Mg II h and k of AB Aur were i nterpreted as 
rotational modulation (e.g. iPraderie et al.l Il986l ICatala et"an 
[l986). Hubrig et al. (2011) detected a rotationally modulated 



* Based on observations collected at the European Southern 
Observatory, La Silla, Chile (Program ID: 084.A-9016, 085.A-9027, 
085.A-9024 



magnetic field of the HAeBe star HD 101412. In this work we 
determine the rotation period of the Herbig star HD 135344B 
by measuring radial velocity variations using multi-epoch high 
resolution optical spectra. 

The paper is organized as follows: Section |2] presents the 
observations and data reduction. The determination of the stellar 
parameters and the radial velocity measurements are presented 
in Sects. [3] and |4] A discussion of the results and conclusions 
can be found in Sects. |5] and |6] The Appendix provides a table 
of measured quantities. 

2. Observations and data reduction 

The observations were carried out in two observing cam- 
paigns w ith the Fiber-fed Ex tended Range Optical Spectrograph 
(FEROS : lKaufer et alJI 19991) at the 2.2m MPG/ESO telescope at 
La Silla Observatory in Chile. FEROS covers the whole optical 
spectral range from 3600 A to 9200 A and provides a spectral 
resolution of *48 000. The fibre apeiture of FEROS is 2 arc- 
sec on sky. A contamination of the spectra due to other stars 
around HD 135344B can thus be ruled out (Sect. |3]l. In total, 
44 spectra of HD 135344B were obtained over a time range of 
five months between March and July 2010. Depending on the 
conditions, we observed the star up to four times a night with 
a separation of two hours between the observations. The aver- 
age exposure time was 16 minutes per spectrum, which results 
in an average signal-to-noise ratio (SNR) of 230 at 5500 A of 
our spectroscopic data set. A set of 40 spectra were obtained us- 
ing the Object-Calibration mode where one of the two fibers is 
positioned on the target star and the other fiber is feed with the 
light of a Th Ar+Ne calibration lamp. This mode allows to mon- 
itor and correct for the intrinsic velocity drift of the instrument. 
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Further four spectra were obtained using the Object-Sky mode 
where the second fiber is pointing to the sky, which allows the 
subtraction of the sky background from the target spectrum. The 
reduction of the raw data was performed using the online data 
reduction pipeline available at the telescop^E The pipeline does 
the bias subtraction, flat-fielding, traces and extracts the single 
echelle orders, applies the wavelength calibration, and corrects 
for the barycentric motion. For each exposure, it produces 39 
individual sub-spectra representing the individual echelle orders 
as well as one merged spectrum. 



3. Astrophysical parameters of HD 135344B 

HD 135344B (SAO 206462) belongs to the Sco OB2-3 
(Upper Centaurus Lupus, UCL) star-f orming region, whos e 
center has a distance of 14 ± 2 pc dde Zeeuw et al.l Il999h . 
iPreibisch & MamaieS (l2008l) also list the individual distances of 
81 group members. The median and standard deviation of these 
values gives 142 + 27 pc for the distance. Note that this spread 
reflects the true extent of Sco OB2-3, rather than observational 
errors in distances to individual stars. We thus adopt a value of 
142 + 27 pc for the distance of HD 135344B in the foUowing. 

HD 135344B is the secondary star of the visual binary 
system of HD 135344 (SAO 206463). The tw o components 
are separated by 21" (PA=197°, .Mason etalj 12001.) . which 
translates to a projected separation of 3000 + 600 AU for the 
given distance of 142 + 27 pc. Therefore, a gravitational inter- 
action between the primary and the disk of HD 135344B can be 
ruled out for fliis large s eparation. Using the H ST/NICMOS2 
camera at A - 1.6 um. lAugereau et aLr(l200ll ) found a close 
binary system lying 5.8" from HD 135344B. A companionship 
was ruled o ut by the detection of different proper motions by 
iGradv et"aD (f2009) which used the same instrument. 

The disk around HD 135344B was subject to numerous 
studies carried out over a large wavele ngth ran^e usinp imaguig, 
spectroscopy and i nterferometry (e.g. Thi et al. 2001; Dent et al. 
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..1,2006: Brown et al. 



2007^,^^ 
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Fedele et al.j,2008^ 



Pontoppidan et alJ 120081: iGradv et all l2009t IVerhoeff et al 
2010) and revealed a complex structure. An outer radius of 
200 AU for the dusty disk was derived by Doucet et al. (2006) 
at 20.5 iJ.m. From modeling the Spitzer spectrum. Brown et al.l 
(I2OO7I) found a gap in the dusty disk between 0.45 AU and 
45 AU. We should note, that lBrown et all (l2007h do not resolve 
the inner disk but estimated the location of the inner rim of 
the gap based on potentially ambiguous SED model fits. This 
is in a greement with the mid-lR observations of IVerhoeff et alJ 
(l2010l) . One explanation for the presence of a gap in the disk, 
besides grain growth and photo-evaporation, is the formation of 
a (su b-)stella r companion which causes a dynamical clearing. 
iPontoppidan et alJ (l2008h rule out the presence of a stellar 
companion around HD 135344B based on their detection of 
molecular CO gas (A = 4.7 yum) at 0.3-15 AU. 
The value for the disk inclination derived from different observa- 



tions at_differentwavel^ va ries from 1 1 ° to 61° dDent et al 
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Fig. 1. A pictographic sketch of the HD 135344B system. The 
star, the dusty disk (brown color) with its large inner gap, the 
gaseous disk (blue), and an hot spot emerging from the inner 
rim of the inner disk caused by active accretion (red) are shown. 



2009h . We should note, however, that differ- 
ent methods probe different spatial scales in the disk. Therefore 
those different values may not be in contradiction to one another 
but may reflect the presence of a more complex morphology 
than assumed in the simple symmetric models by which these 



values are derived. iGradv et"al] (1200 9') ruled out inclinations 
greater than 20° for the outer disk which is in good ag r eemen t 
with the value of 11° + 2° measured by Ibent et al I (l2005h . 
Therefore, the outer disk is se en almost face-on as it is already 
indicated by direct imaging (iDoucet et al.l l2006t iGradv et"aD 
2009). Interferometric N-band observations of the inner part 
of the star-disk system revealed the presence of a much higher 
inclined (53°-61°) structure with an estimated in ner radius of 
about 0.05 AU and an outer radius of 1.8 AU (iFedele et al.l 
2008). 

The disk of H P 135344B con tains (2.8 + 1.3) ■ 10"^ Mq of 
gas and dust dThi et al.l 12001 1) . The dete ction of polycyclic 
arom a tic hydrocarbon (PAH) features (e.g. ICoulson & Walthed 
1 19951: iBrown et all 120071: iFedele et al.l l2008l) implies a rich 
chemis try in the vicinity o f the star. From J=3-2 '^CO measure- 
ments. Dent et aTl (2005) derived an inner radius of <10 AU 
and 75 ± 5 AU for the outer radius where emitting gas is present. 

HD 135344B also shows si gns of ac tive mass accretion. 
This matches the findings of iPott et all (2010), who studied 
young transitional disk systems with a gap in the circumstellar 
dust distribution, similar to HD 135344B (Fig. [TJ. They sys- 
tematically f ound sub-AU scale dust i n accreting transitional 
disk systems. iGarcia Lopez et all d200 6') derived a value for the 
mass accretion rate of Mace ~ 5.4 ■ 10"^ Mp yr'' usin g the Bry 
emission line in the NIR. From FUV data, iGradv et al.l (120091) 
concluded that HD 135344B drives no jet, but found indications 
for a stellar wind leading to a mass-loss rate between 10"" to 
10-**Mayr-i. 

Summarizing all results from the previous observations, 
the following picture of the disk geometry can be drawn: dust 
is present within an inner ring from 0.05 AU to 1.8 AU. This 
structure is highly inclined compared to the outer disk. There 
are indications that mass accretion from the inner disk onto the 
star is still present. Between 1.8 AU and 45 AU, a large gap 
is present. At 45 AU the outer disk begins and is extended to 
about 200 AU and is seen almost pole-on. See Figure [T] for a 
pictographic sketch. 

3.1. Stellar parameters 



Table[T]lists the derived stellar parameters of HD 135344B. The 
luminosity, L*, was computed from integrating the total flux un- 
der the SED. Teff, \ogg, and vsin/* were computed by a self- 
developed tool for retrieving stellar parameters of Herbig stars 
based on fitting synthetic spectra to the observed stellar spec- 
' |http : //www ■ eso . org/sci/f acilities/lasilla/instruments/tfamsT he computation of the synthetic spectrum was carried 
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Fig. 2. The spectral windows used for fitting the synthetic spec- 
trum (red line) to the observed stellar spectrum (black line, 
JD=2455347.63521)ofHD 135344B. 

Table 1. Stellar parameters of HD 135344B. 



Parameter 




Value 


distance 




142 ± 27 pc 


iogL^, 






Tts 




6810 ± 80 K 


iogg 


4.4 


± 0.1 cm s~^ 






1.7t? Me 


age 




9 ± 2 Myr 




1 


.4 ± 0.25 Ro 


<RV^ > 


2.5 


± 1.5 km s"' 


V sin i-i. 


82.5 


± 2.9kms-' 



Notes. < RV^, > represents the mean value of the observed RVs. The 
error does not present the accuracy of the determined /?V* but reflects 
the scatter on RV^, around the mean value. 



out using sPECTRUivQ (iGrav & Corballv 1 19941) with the atlas9 
atmosphere models (ICastelli & Kuruca 2004). Fo r a detailed 
description of this procedure we refer to Miiller et al.l (1201 ll 
in preparation). Figure |2] shows a part of an stellar spectrum 
of HD 135344B (black line, observed at JD=2455347.63521) 
which is well fitted by a single star synthetic spectrum (red line). 
From the position of HD 135344B in the H-R diagram (Fig. [3]), 
we obtained the stellar mass and radius and its age. The derived 
parameters are vy ell in a greement wit h htera ture values, e.g., 
Ivan der Plas et al.' feOOS): Manoi et al.l (l2006l) : [van Boekel et all 
(2005J; Dunkinetal. (1997). The measured effective tempera- 
ture is comparable to that of an F3Ve to F4Ve main-sequ ence 
star and agrees with the value found bv lDunkin et al.l (Il997h . 
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Fig. 3. The position of HD 135 344B in an H-R d iagram along 
with the evolutionary tracks of ISiess et al. I (l2000l) for a metal 
abundance of Z=0.02. 



makes a determination of RVi, difficult since the accuracy is lim- 
ited to some 100 m s"'. For the cross-correlation, only the spec- 
tral range between 4000 A to 7875 A was considered because 
the instrument efficiency drops significantly outside this spec- 
tral range. In addition, areas with strong emission lines, telluric 
lines, and Balmer lines were carefully marked out. The resulting 
cross-correlation function was fitted by a Gaussian function. The 
position of the center of the Gaussian yields RV*. After the in- 
dividual RVi, values were derived separately for each echelle or- 
der, the median value and standard deviation of RV* were com- 
puted. We applied to the RV* values a 1 - cr chpping to remov e 
anomalous velocity values, similar to e.g. iBames et al.l (120051); 
iJefferv et al.l (l2007h . The deficient RV-i, estimates, rejected by 
this cr clipping, do not follow a simple Gaussian noise statis- 
tics but are far off the median RV*, which indicates individual, 
systematic biases of the RV^, estimate in the rejected orders. The 
final RV-i, were then computed based on averaging ten to twelve 
remaining individual orders. The error of the final RVi, is the 
standard deviation of the mean. All measured RV* values are 
listed in Table lA.ll provided in the appendix. Figure |4] shows 
the measured RVs for all 44 spectra with a mean RV* value of 
2.5 km s ', indicated by the gray horizontal dashed line. Peak- 
to-peak variations of up to 2.9 km s ' are present. 



4. Data analysis and results 

4.1. Radial velocity measurements 

The measurement of the stellar RV has been done by cross- 
correlating the stellar spectrum with a template spectrum. The 
template was a synthetic spectrum representing the fitted stellar 
parameters of HD 135344B (Table [TJ. We measure a projected 
rotational velocity, v sin of 82.5 km s"' (Sect. lTTT ). Therefore, 
all stellar spectral lines are highly broadened and blended, which 

^ |http : //www ■ phys . appstate . edu/ spectrum/spectrum . html | 



4.2. Radial velocity variations 

In order to identify periodicities present in the RV-t data, we 
computed the generalized L omb-Scargle (GLS) periodogram 
(IZechmeister & Kursterll20()9) and its window function (Fig.|5]l. 
The GLS exhibits a strong peak at a period of 0. 16045 d (in the 
following we use 0.16 d), marked by the red arrow. Its false- 
alarm probability (FAP) is 0.64%. In addition, there are several 
other peaks in the GLS present which are less or not significant 
compared to the 0.16 d period. This makes the interpretation of 
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Table 2. Results of the sinusoidal fit to the RVi, data. 
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Fig. 4. Measured RVs for HD 135344B. The horizontal gray 
dashed line represents the mean value of the measurements. 



a periodogram ambiguous. In order to identify alias frequencies 
caused by uneven sampled data in time, the window function 
has to be considered. Significant peaks in the window function 
can produce aliases in the GLS which are separated from the 
true peak by the frequency difference. Figure|6]shows a close-up 
view of the GLS, centered around the 0. 16 d period (or at the fre- 
quency of 6.232 d"' respectively). The two strongest other, but 
less significant peaks are located at 5.230 d"' and 7.234 d"' . The 
difference of both frequencies to 6.232 d"' is 1.002 d"', which 
corresponds to the sidereal day. Figure |7]b) shows a close-up 
view of the window function centered around the significant fre- 
quency 1 .0 d"' . The highest peak is at a frequency of 1 .002 d ' . 
Therefore, we can conclude that the peaks to the left and to the 
right of the 0. 16 d period in the GLS are aliases. Each of the three 
peaks in Fig.|6]show side lobes which are separated by 0.021 d"' 
for all cases. From the close-up view of the window function 
(Fig. I?] a), we can identify a significant frequency present at 
0.021 d ' which is exactly the measured frequency difference 
between the side lobes of three main peaks. The 0.021 d ' fre- 
quency corresponds to a period of about 48 d. By looking at the 
distribution of observations in Fig.|4] we can identify three sep- 
arated data sets (the first data set consists of the first single data 
point only). The time difference between the first and the second 
data set is about 100 d, and the difference between the second 
and the third data set is about 50 d. This explains the observed 
aliases at this frequency difference in the GLS and the presence 
of a significant peak in the window function. The observations 
cover a time range of about 150 d. Therefore, the window func- 
tion (Fig.|2]a) shows a significant peak for frequencies smaller 
than 0.006 d"' (or periods greater than 150 d respectively), but 
they are not visible as side lobes of the main peaks in the GLS. 
In addition, we subtracted a sinusoidal fit from the RV* data and 
computed a GLS periodogram of the residuals (Fig. |9) in order 
to verify that the 0.16 d period is present in the RVi, data. The 
GLS periodogram of the residuals clearly shows that the peak 
at 0.16 d as well as its corresponding aliases disappeared. There 
are no further significant peaks at other periods present. From 
this analysis, we consider the 0. 16 d period as significant and as 
the only real period present in the observed RVi, data. 
Figure |8] shows the phase-folded RVi, data for the 0.16 d period. 
The solid line represents a sinusoidal fit to the data with a period 
of 0.16045 + 2 ■ 10"^ d and an amplitude of 493 ± 86 m s"'. The 
results of the fit are presented in Table |2] 



Parameter 



Value 



rms 
period 
amplitude 
phase (to JD„,in) 



2.5 

478 m S-' 
0.16045 ±2 - 10^^ 
493 ± 86 m s"' 

0.95 ± 0.03 
2508 ± 61 m s"' 
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Fig. 5. GLS periodogram (upper plot) and window function 
(lower plot) of the RVi, data. The significant period at 0.16 d is 
marked by the red arrow. FAP thresholds for 0. 1 %, 1 %, and 10% 
are indicated by the horizontal dashed lines. The data are plot- 
ted over frequency but the upper x-axis gives the corresponding 
period. Note that periods greater than 1 d are in the frequency 
range between and 1 d 



4.3. Bisector analysis 

Stellar activity like cod and hot spots, granulation and pulsation 
can mimic significant RV variations that can be periodic due to 
stellar rotation. In case of a cold (i.e. dark) spot, the line profile 
becomes asymmetric because of the reduced stellar light. Stellar 
rotation leads to a modulation of the asymmetry and causes RV 
variations depending on the stellar inclination and latitude of the 
spot. The analysis of the line profile (in the following bisector) 
is routinely applied in the search for extrasolar planets when the 
RV technique is used in order to identify if the observed pe- 
riodic RV variations are caus ed by stellar activity or due to a 
companion (e.g. Oueloz et al. 2001). A correlation between RV^, 
and bisector velocity span would indicate the the RV* variation 
is caused by line shape asymmetries, i.e. by rotational modu- 
lation, e.g., due to a spot. The shape of the computed cross- 
correlation function (CCF) represents the mean line profile of 
the selected lines in the observed spectra. We measured the bi- 
sector velocity span (BVS), bisector displacement (BVD), and 
bisecto r curvature (BC) fo r all our spectra following the defini- 
tion bv'Povic h et all (HooT). To construct the bisector, the middle 
points of the horizontal segments, which connect a point on the 
left and on the right side of the CCF at the same flux level, are 
computed. The bisector gets divided into three zones. The BVS 
is simply the difference between the upper and lower zone. The 
BC is defined as the difference between the upper and lower bi- 
sector spans. To derive the BVD the average of the three bisector 
velocity zones is computed. For the computation of the individ- 
ual quantities we chose heights of 30%, 50%, and 85% of the 
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Fig. 6. Close-up view of the GLS centered around the most sig- 
nificant peak at 0.16 d. Two less significant peaks at 5.230 d"' 
and 7.234 d"' are present in addition and are aliases caused by 
the sidereal day. All three peaks show side lobes with a differ- 
ence of 0.021 d"' in frequency, caused by the uneven sampling 
of our data. 
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Fig. 7. Close-up views of the window function for the significant 
peaks at <0.006 d"', 0.021 d"', and 1.002 d"'. 



bisector Similar to the RVi, measurements, we measured BVS, 
BVD, and BC for each computed CCF. The final bisector val- 
ues were derived by averaging the single values. The errors of 
BVS, BVD, and BC are the standard deviations of the means. 
In Figure [To] we plotted all three bisector quantities against RV. 
The computed linear correlation coefficients r are noted in the 
upper right corner of each plot. For none of the three quantities a 
significant linear correlation is present. Only a weak trend might 
be present for BVS and BD. However, the sensitivity of the bi- 
sector method d ecreases rapidly w hen it comes to low stellar 
inclinations (e.g. lDesort et al.ll2007h . For HD 135344B it can be 
assumed that we see the star almost pole-on (Sect. |5]l and there- 
fore we would not expect a strong correlation between BVS and 
RV. 




0.0 0.5 1.0 1.5 2.0 

Phase 



Fig. 8. Phase-folded plot of the RVi, data. The solid line is a si- 
nusoidal fit with a period of 0. 16045 ±2 ■ 10"^ d and an amplitude 
of493 + 86ms-i. 
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Fig. 9. GLS periodogram of the residual RV* data after subtrac- 
tion of the sinusoidal fit. The red arrow marks the position of the 
0.16 d period which is no longer present, as well as its corre- 
sponding aliases. The new peak with the highest power at 0.14 d 
has a FAP of 3.4% and is therefore not significant. 

4.4. Ha measurements 

The Ha line profi le of HD 135344B und ergoes rapid and sig- 
nificant changes. iFinkenzeller & Mundtl (11984.) classified the 
Ha line profile of HAeBe stars into three categories: single- 
peak, double-peak, and P Cygni profiles. The Ha line profile 
of HD 135344B adopt each of these types at different times. 
Figure [TT] shows the Ha line profiles of our spectroscopic data 
set of HD 135344B ordered with respect to their phase value 
according to a 5.77 d period. A double-peaked profile indi- 
cates accretion. A blue-shifted absorption (P Cygni profile) indi- 
cates an outflow of material. We measured the equivalent width, 
^^(Ha), of HD 135344B in ourFEROS spectra. The measured 
values are listed in Table [ATI provided in the appendix. For the 
measurement of ^^(Ha), we selected several distinct wave- 
length intervals on both sides of the Ha line and used a 2nd or- 
der polynom to describe the continuum. This was repeated four 
times for different intervals. The deviation of all four measure- 
ments defines our error For ^^(Ha), we obtained a mean value 
of -10 A. The GLS periodogram of the £'W(Ha) values did not 
show a significant period. 

We also measured the full width of Ha at 10% height, 
Wio(Ha), which is frequently used for TTSs in order to dis- 
tinguish between chromospheric activity and accretion of cir- 
cumstellar material CWh ite & Basri 2003; Natta et al.l l2004l 
iJavawardhana et al.l2006l) . We applied the same method for nor- 
malizing the spectra as for EW(lia). We measured a mean value 
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Fig. 10. Bisector quantities vs. RV-^. For each plot the linear cor- 
relation coefficient r and the proba bility p that the data are lin - 
early uncorrected are computed (IBevington & Robinsonll2003l) . 
There is no linear correlation between the bisector quantities and 
RVi, present. 



of 9.8 A for H^ioCHo-), which correspond s to an accretion rate of 
3 ■ 10 '' Mo yr ' if Eq. 5 in lNattaet al.1 (l2004 is appHed. The 
measured values of Wi()(T[q;) are listed in Table I A. fl provided in 
the appendix. The periodogram of Wio(Ha) shows a significant 
period at 5.77 d with a FAP of 1.9 ■ 10"^. A phase-folded plot for 
the Wio(Hq') data (black data points) and the sinusoidal fit (red 
line) is shown in Fig. [12] Assuming Keplerian rotation of the 
circumstellar gas, this period corresponds to a radial distance of 
^ 0.07 AU. The error bars of Wio(Hq') are the statistical errors of 
the single measurements and reflect the precision of the individ- 
ual measurements. These error bars are most likely overestimat- 
ing the accuracy of each measurement, since the distance of the 
data points to the sinusoidal fit often exceeds the error bars sig- 
nificantly. A possible reason for an accuracy bias of individual 
Wio(Hq') measurements can be that the Ha line can form in dif- 
ferent locations. Also, the formation can be caused by different 
processes, e.g. accretion column, accretion shock on the stellar 
surface, disk surface, and wind. Therefore, we see a significant 
scatter of the measured Wio(Ha') values around the sinusoidal 
fit, the level of which exceeds the precision of the individual 
measurements. However, the bias of each individual measure- 
ment statistically spread around the fit. We binned the phase in 
0. 1 wide bins and averaged the Wi()(Hq') values lying in the cor- 
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Fig. 11. Ha line profiles of HD 135344B. The plots are ordered 
with respect to their phase value according to a 5.77 d period, 
which is shown in the upper left comer of each plot. 



responding bin. The shown error is the standard deviation of the 
values in each bin. These values are overplotted in Fig. [T2]with 
green triangles. The green error bars are now reflecting both the 
precision and accuracy of the Wio(Hq;) measurements and are 
well described by the sinusoidal fit. 

5. Discussion 

The observed 0. 16 d period might be caused by rotational modu- 
lation of one or more spot(s) on the stellar surface. Adopting the 
stellar radius of 1 .4 + 0.25 R© and v sin /* = 82.5 + 2.9 km s"' 
from Table [1] the inclination of the stellar rotation axis is = 
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11.0 




Fig. 12. Phase-folded plot of the Wio(HQf) data. The solid line is 
a sinusoidal fit with a period of 5.77 + 0.02 d and an amplitude of 
0.70 + 0.08 A. The green triangles are the average values com- 
puted by binning the phase in 0. 1 wide bins. Their error bars are 
derived by computing the standard deviation of the correspond- 
ing Wio(Ha) values in each bin. 



11° ±2°. We thus see the star almost pole-on. This could explain 
the weak correlation between the BVS and RVi, (Sect. 1431 ). The 
disk of HD 135344B is also foun d to be almost face -on with in- 
ch nation meas urements of < 20° dGradv et al.l2009h and 11° +2° 
bv lDent et all (1200 5 ). i.e., the orientation of the stellar rotation 
axis, derived from the assumption that the 0. 16 d RVi, period rep- 
resents the stellar rotation period, would be almost perpendicular 
to the mid-plane of the outer disk. This supports the hypothesis 
that the 0.16 d period is caused by photospheric effects. 
Another approach to verify the origin of the 0.16 d period is to 
apply a Keplerian fit to the data. Under the assumption of a cir- 
cular orbit (eccentricity is set to 0) and a period of 0.16 d, the 
semi-major axis is 7-10"'' AU, which corresponds to 1.5 R©. For 
HD 135344B, we found the stellar radius to be 1.4 + 0.25 Rq 
(Table [1]), i.e, the semi-major axis would be in the range of the 
stellar radius. This shows that the companion hypothesis is not 
plausible, and thus favors a photospheric interpretation. Thus, 
the rotational period of 0.16 d or 3.9 hours makes HD 135344B 
an extremely rapid rotator 

5.1. HD 135344B is rotating close to break-up velocity 



The break-up velocity Vc - yfUMZlR* of the star, where G is 
the gravitational constant, has a value of 480 + 60 km s"^ Using 
the value of = 1 1° + 2° derived from the stellar parameters 
leads to a true rotational velocity of 432 + 81 km s"' . Therefore, 
HD 135344B is rotating at or close to the break-up velocity at its 
equator. 

[Aufdenb erg et al.l (l2006l) measured the polar and equatorial ra- 
dius of Vega interferometrically which rotates at 91% of its 
break-up velocity. They measured a difference of 20% between 
equatorial and polar radius and found the temperature gradient 
to be 2250 K. For HD 135344B, similar effects are expected. 
Because the star is seen almost pole-on, we might overestimate 
its luminosity and thus its mass and temperature. The angular 
stellar diameter of HD 135344B is about 100 /ias and not resolv- 
able by current stellar interferometric facilities which prevents a 
direct confirmation and the determination of its true stellar pa- 
rameters. 

The continuous accretion of circumstellar material onto the star 
could have caused the spin-up of the stellar rotation up to v^.. 



HD 1 35344B has only a very weak magnetic field (iHubrig et alJ 
l2009l) . Therefore, the star is probably decoupled from its disk, 
which prevents the star from magnetic braking and loosing spin. 
In contrast, for the majority of TTSs only moderate vsin/* val- 
ues are observed which indi cates a disk-brakin g mechanism dur- 
ing the accretion phase (e.g. lWeise et al.ll2010h . 

5.2. Possible origins of the observed periodicity 
5.2.1. The von Zeipel effect 

Because of its rapid rotation, the star should show an oblate 
shape of its stellar atmosphere due to the centrifugal force. This 
results in a redistribution of the flux whic h is proportional to 
the local surface gravity (Ivon Zeipellll924l) . i.e., the equatorial 
zones are darker with a lower Te/f than the polar zones. This 
phenomenon is known as gravity-darkening and is symmetric 
with respect to the equator In order to produce a periodic RV^, 
signal, i.e. mimic a stellar spot, the flux and temperature gradi- 
ent has to be off-centered from the rotational axis in order to pro- 
duce a Doppler shift, i.e. the von Zeipel effect cannot account for 
the observe d periodicity of the spectra. Smith & Wo rlevI (Il974l) 
and ICIare^ ([2000) showed that for differentially rotating radia- 
tive stellar atmospheres, the flux over the stellar surface varies 
stronger than the surface gravity, which could produce asym- 
metric flux distributions. 



5.2.2. Stellar spots 

We computed different spot configurations (temperature differ- 
ence of the spot with respect to the photosphere Arspot, filling 
factor /spot, latitude of the spot 0) taking the determined stellar 
parameters into account. 

The effects of star spots on photometry and RV* were numer- 
ically derived with a relatively simple model. A sphere with 
radius R, unit surface brightness /, and one or more round spots 
with spot filling factor /spot at specified latitude and longitude 
was projected on a 2D Cartesian grid with 200 x 200 pixels. 
The relative spot brightness at the considered wavelength was 
calculated from the ratio of effective photospheric temperature 
to spot temperature, assuming black-body radiation. To account 
for opacity effects on the stellar surfa ce, a linear limb darkening 
law with coefficients adopted from Claret & Hauschildi (l2003h 
was applied. The inclination of the rotation axis can be arbi- 
trarily chosen. Solid-body rotation was assumed; differential 
rotation was not considered. Each cell of the star image was thus 
assigned a brightness and a radial velocity value. The surface 
brightness distribution at each phase of the stellar rotation (in 
steps of 1 deg) was then integrated to obtain the total brightness 
and the position of the photo center (P' moment). The P' 
moment of the radial velocity distribution, weighted with the 
surface brightness distribution, was adopted as a proxy for the 
the mean radial velocity. Absorption line shapes were not taken 
into account. 

The two extreme scenarios which can explain the RVi, variations 
are: 1.) dark spot, Arspot = -1350 K, = 40°, /pot = 5%, 2.) 
bright spot, AFspot = 2000 K, = 85°, /pot = 56%. The large 
covering fraction of the latter model may more closely resemble 
the von Zeipel effect than traditional star spots. So far, there 
are no photometric monitoring data for HD 135344B available 
which would set additional constraints on the possible spot 
configurations. 

We conclude that both explanations for the origins of the 
observed RV-t variations, such as stellar spots or the flux 
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and temperature gradient between the pole and equator, can 
adequately explain the data. 

Due to the accretion process, the Wio(Ha) variations may 
trace a hot spot from an accretion funnel flow which has its 
origin at the inner edge of the inner disk. A periodic variation of 
5.77 d is found for Wio(Hq') (Sect. l4~4l ). This period corresponds 
to a Keplerian radius of a!0.07 AU th at is close to the inner edge 
of the disk estimated to 0.05 AU bv lFedele etalj (|2008), given 
the uncertainty. A similar effect w as observed for EW (lia) 
variations for the TW Hya system bv lSetiawan et al.l (l2008h . 
Several studies have shown that the disk lifetim es of young 
stars are expected to b e ~5 to ~10 Myr (e.g. | Haischetal] 



20011: iB ouwma n et alJ [20061: Ijavawardhana et al.1 l2006t 



Siciha-Aguilaretal.1 l2006t lFedeleet al.1 I2OI0I ). At this age. 



stars lose their inner dusty disk and accretion stops. This seems 
to coincide with t he dissipation of gas in the inner disk, too 
dJavawardhana et al., .2006) . In contrast to that, HD 135344B 
shows still accretion signatures, an inner dusty disk with an 
inclination, which is much greater than that of the outer disk, as 
well as a gas-rich inner environment at an age of 9 ± 2 Myr A 
possible explanation might be that HD 135344B is losing parts 
of its outer atmosphere due to the rapid rotation at or close to its 
break-up velocity and feeding the inner part of the disk with its 
own material, which then accretes back onto the star. 
A possible stellar or sub-stellar compa nion located in th e gap of 
the disk suggested by literature (e.g. iGradv et al.ll2009b cannot 
be confirmed or ruled out with our measurements because our 
spectroscopic data set covers only five months in total and 
the uncertainties in the RVi, measurements are of the order of 
300 ms-'. 



6. Conclusions 

In this paper we presented the first direct measurement of the 
stellar rotational period of a Herbig star The period was de- 
rived from RVi, measurements by cross-correlating optical high- 
resolution spectra with a synthetic spectrum. The rotational pe- 
riod of HD 135344B was found to be 0.16 d. In addition, we 
determined new reliable stellar parameters for HD 135344B, 
independent from previous literature values, and demonstrated 
the possibility to measure the RVi, of a Herbig star with a high 
V sin ii, value. From the measured stellar parameters, we were 
able to estimate the stellar inclination to a value of 11° ±2°. 
With these data, we concluded that HD 135344B is rotating at or 
close to its break-up velocity and feeding its stellar vicinity with 
gas and dust. This is the first observational evidence that a young 
intermediate-mass star can rotate close to its break-up velocity. 
The direct determination of the stellar rotational period by RVi, 
measurements may be applicable also to other Herbig Ae/Be 
stars. A larger sample could therefore reveal the evolution of 
angular momentum, the presence of star-disk coupling mecha- 
nisms and the presence possible magnetic fields of these young 
intermediate-mass stars. 

We also found variations in the IVio(Ha') measurements with a 
period of 5.77 d, corresponding to a Keplerian radius close to the 
inner edge of the inner disk of HD 135344B, which could be a 
hot spot caused by an accretion funnel flow. 
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Table A.l. Measured /JV*, EWiYla), WioCHor), and their corresponding uncertainties. 



Time 

[JD-2400000d] 


[m s-'] 


0■RV^, 

[m s-'] 


EWiHa) 
[A] 


C EW(Ha) 

[A] 


Wio(Ha) 
[A] 


O" IV(Ha) 

[A] 


55256.86690 


2645 


335 


-6.02 


0.21 


9.80 


0.07 


55346.70410 


1593 


302 


-9.84 


0.13 


9.14 


0.05 


55347.63521 


1859 


315 


-10.44 


0.17 


8.54 


0.06 


55348.73788 


2041 


260 


-13.17 


0.19 


9.24 


0.02 


55349.74833 


2223 


400 


-11.42 


0.20 


8.45 


0.08 


55351.69950 


2331 


112 


-12.80 


0.24 


8.45 


0.03 


55355.49994 


2108 


192 


-8.83 


0.12 


8.37 


0.08 


55355.68805 


2300 


195 


-8.96 


0.13 


8.73 


0.10 


55355.85300 


1784 


784 


-9.85 


0.20 


8.84 


0.11 


55356.78750 


1529 


303 


-8.24 


0.20 


9.35 


0.07 


55358.63746 


3925 


494 


-9.17 


0.15 


9.20 


0.03 


55358.71724 


2364 


408 


-9.15 


0.14 


9.21 


0.06 


55360.72711 


2664 


358 


-10.35 


0.13 


8.86 


0.13 


55362.48448 


2432 


215 


-9.55 


0.15 


7.78 


0.04 


55362.57222 


1732 


205 


-9.67 


0.11 


7.97 


0.07 


55396.51686 


3415 


245 


-10.48 


0.20 


7.63 


0.06 


55396.57497 


2641 


306 


-10.73 


0.08 


7.90 


0.03 


55396.61357 


2388 


123 


-10.91 


0.17 


8.37 


0.02 


55396.66325 


3566 


214 


-10.90 


0.16 


8.69 


0.05 


55396.69535 


2672 


382 


-11.18 


0.08 


8.75 


0.03 


55397.52834 


1992 


339 


-10.01 


0.11 


8.27 


0.05 


55397.59554 


2288 


247 


-9.99 


0.13 


8.11 


0.11 


55397.64898 


2739 


381 


-10.02 


0.30 


8.05 


0.03 


55397.68781 


2866 


241 


-10.27 


0.15 


8.06 


0.09 


55398.51102 


1932 


519 


-11.25 


0.29 


8.45 


0.09 


55400.51882 


2692 


303 


-10.66 


0.14 


10.03 


0.06 


55400.57368 


2159 


173 


-10.18 


0.22 


9.65 


0.06 


55400.61784 


2891 


239 


-9.70 


0.14 


9.64 


0.10 


55400.67505 


2803 


321 


-9.52 


0.20 


9.41 


0.07 


55401.48617 


3341 


376 


-10.65 


0.17 


9.20 


0.03 


55401.63979 


3461 


215 


-10.01 


0.17 


9.23 


0.02 


55401.69736 


1035 


546 


-9.56 


0.21 


8.93 


0.07 


55402.56980 


2032 


374 


-9.41 


0.00 


8.63 


0.06 


55402.62734 


3849 


347 


-9.35 


0.24 


8.36 


0.11 


55402.69271 


1769 


433 


-8.77 


0.18 


8.09 


0.06 


55403.56511 


2802 


378 


-8.49 


0.19 


8.81 


0.04 


55403.59827 


2209 


193 


-8.08 


0.23 


8.48 


0.13 


55403.63910 


2810 


471 


-8.05 


0.19 


8.60 


0.11 


55403.67692 


1898 


311 


-7.99 


0.23 


8.59 


0.11 


55404.50509 


3058 


250 


-8.80 


0.18 


9.55 


0.05 


55404.67289 


2932 


247 


-9.56 


0.18 


9.61 


0.06 


55405.63121 


2431 


197 


-12.65 


0.16 


9.74 


0.07 


55405.68095 


2686 


313 


-12.62 


0.19 


9.55 


0.04 


55407.58835 


3106 


371 


-10.67 


0.16 


8.52 


0.02 
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